Exotic pairing via a central attraction: from the BCS to the Bose limits 
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In the context of a simple model featuring an explicit, central interaction potential, and using a 
standard functional-integral technique, we study superconductivity with angular momentum quan- 
tum number I = 2 as an emergent property of the many-body system. Our interaction potential 
is attractive at a finite distance ro, and the breaking of the rotational symmetry is the result of 
an interplay between ro and the inter-particle distance r s which we deem generic to interactions of 
this type. However such interplay, responsible for the preference of a d-wave state for a range of 
intermediate densities, takes place only in the BCS limit. In contrast, as the Bose-Einstein (BE) 
limit is approached the internal energy of the "preformed pairs" becomes the dominant contribution 
and there is a quantum phase transition in which the s-wave symmetry is restored. We also find 
that the limiting value of the critical temperature is fcsT c — ► 3.315 h 2 /2m* [n/2 (21 + l)] 2 ^ 3 , which 
coincides with the usual result only for I — 0; for I > 0, it differs in the degeneracy factor 1/ (21 + 1), 
which lowers T c . Our results thus place constraints on exotic pairing in the BE limit, while at the 
same time indicating a particularly interesting route to pairing with I > in a BCS superconductor. 



I. INTRODUCTION 

It is a surprising consequence of BCS theoryi that 
certain central interaction potentials V (|r — r'|) lead to 
Cooper pairing with a finite value of the angular mo- 
mentum, thus breaking the rotational symmetry of the 
continuumH. The phenomenon is analogous to what 
happens in "Hubbard" models with attraction between 
nearest-neighbours for which, as is well knownQ, d- 
wave pairing can break the symmetry of the crystal lat- 
tice. However, it is especially interesting in the original 
contexta of a continuum modpl, because of the contrast 
with the well-known theoremta for two-body pairing in 
real space, which demands that the ground state min- 
imizes the orbital angular momentum. This rotational 
symmetry breaking is thus a many-body effect, and one 
expects that in the limit of low densities and strong at- 
traction, when the BCS ground state is a jBoaft-Einstein 
(BE) condensate of non-overlapping pairsBErQIa, the ro- 
tational symmetry of the system is restored. However 
until now investigations of the BCS to Bose crossover 
for non s-wave pairing were performed either in the 
context of lattice models (see Ref. |^ for a review and 
Refs. ^ 10 n],|l^|L3 14 for some examples of recent wopk) 
or for the anisotropic interaction potential of Ref. [l5|.Efl 

In this paper we take a slightly different approach 
by studying a continuum model, but choosing to work 
with an explicit, central interaction potential V (|r — r'|) 
which can lead to pairing with more than one vahie.of the 
angular momentum: the "delta-shell" potentiaOO. The 
resulting "delta-shell" model (DSM) captures, in an ide- 
alised way, the essential feature leading to Cooper pairing 
with finite value of the angular momentum, namely be- 
ing attractive at a well-defined, finite distanceEJ. Thus 
we expect some of the novel features that we shall de- 
scribe, pertaining to the mechanisms by which the rota- 
tional symmetry is broken in the BCS limit and restored 
in the BE limit, to be generic to a large class of cen- 



tral effective interactions. In particular we shall see that 
for such models the evolution of a BCS superconductor 
with exotic pairing towards the BE limit involves a phase 
transition in which the symmetry of the superconducting 
order parameter is increased. This adds to the work by 
Babaev and KleinerttSl who also found, in the context of a 
chiral Gross-Neveu model, a phase transition associated 
with the BCS to Bose crossover. However the nature of 
the phase transition that we describe here is quite differ- 
ent, as it takes place in the superconducting state, while 
that of Babev and Kleinert corresponds to the formation 
of preformed pairs in the normal state. 



II. THE DELTA-SHELL MODEL 

The first discussions of exotic Cooper pairingi took 
place in the context of the weak-coupling theory of su- 
perfluid 3 He. It was assumed that there existed a central, 
non-retarded interaction potential V (|r — r'|) acting be- 
tween particles at positions r and r'. One then writes 

oo 

V (k - k') = J2 K v d k l . l k 'D (21' + 1) Pv (k-k') , (1) 

l'=0 

where 7 (k-k') = / d 3 r e i ( k - k ')- r V (r - r'), and finds 
that each of the terms in this series leads to pairing with a 
different value of the angular momentum quantum num- 
ber, I. As it can, and has been, argued, in the weak- 
coupling limit one can approximate 



V (k-k') ^Ki (21 + 1) Pt (k.k' 



(2) 



where I is the value of I' for which the coupling constant 
on the Fermi surface, 



K v = K v (k F , k F ) , 



(3) 



2 



is largest. The approximate form (J2J) of the potential 
V (k — k') is, for I > 0, anisotropic, and it leads to pair- 
ing with finite angular momentum quantum number I 
(see Ref. 0). For I = 0, it reduces to the BCS "contact 
potential"!!, leading to s-wave pairing. Although intro- 
duced in the context of a weak-coupling theory, the con- 
tact potential h^^^ieED-Josed-tQ study the BCS to 
Bose crossoverBBMBHeaO-BBo. Similarly, Stintz- 
ing and Zwerger have considered a simplified potential 
of the form || with I = 2 to study the^CS to Bose 
crossover for pairs with d x 2_ y 2 symmetry!^ (but in two 
dimensions, and with the additional assumption of sepa- 
rability to make itpHiore tractable) . One of the key results 
of this later workta was that the critical temperature is 
given, in thepBE, limit, by the same expression as in the 
s-wave caseEHHirE3: 

h 2 /n\ 2 / 3 

UbT c ~ 3.315 I — ) for s andd x 2_ y 2 pairing (4) 

Although very useful, the above approach is not ap- 
propriate to study the question that we are interested in 
here, since it introduces a particular pairing symmetry 
at the level of the interaction potential. In contradistinc- 
tion, we want to find pairing with I > as an emergent 
property of the many-body system. Moreover, we would 
expect, on the basis of the above arguments, to recover 
/ = pairing in the BE limit, in which the internal struc- 
tures of the Cooper pairs are independent. This physics 
seems also to be absent from those studies, as Eq. (Q) 
suggests that the critical temperature is degenerate for s 
and d x 2_ y 2 superconductivity. 

An alternative strategy is to do the calculations tak- 
ing the full r-dependence of V (r) into account. A-study 
of this type was carried out by Andrenacci et a/Jlil who 
took a Gaussian form for V (r) . This allowed them to 
investigate the properties of the crossover at finite den- 
sities (in contrast, as is well known, the procedure re- 
quired to regulate the ultraviolet divergences associated 
with simplified potentials of theJorm Q at all couplings 
is only valid in the dilute limitEa). They also considered 
the highly idealised separable potential introduced in the 
seminal paper by Nozieres and Schmitt Ifew, which has 
been employed in several other instancesc3'E3 on account 
of its mathematical simplicity (but note that this is not, 
strictly speaking, a central potential). However in either 
case there was no rotational symmetry breaking: even at 
high densities they only obtained s-wave superconductiv- 
ity (the discussion of d x 2_ y 2 superconductivity in Ref. [ll] 
is based on a lattice model). _ . 

On the other hand, a simple argument!!^ based on the 
BCS "gap equation" suggests that I > Cooper pairing 
is associated with central potentials V (r) that are non- 
monotonic functions of r, with maximum attraction near 
some finite distance, r ~ ro > 0. The delta-shell poten- 
tial was proposed in Ref. [ll] as the simplest form of V (r) 
that has this feature: 

y(|r-r'|)=-^(|r-r'|-ro) (5) 




r 



Figure 1: The delta-shell interaction potential. Left: the two 
particles attract each other only when each of them lay on 
a thin shell, of radius ro, centered on the other one. Right: 
The delta-shell interaction potential can be regarded as an 
approximation to any central potential that is attractive only 
near some distance ro (see text). 

The resulting DSM can be regarded as the continuum 
analogue of the lattice model with nearest-neighbour at- 
traction discussed in Refs. |^, |Io| , pT|Jl^Jl^ , for example. 
But note that in the DSM the distance ro at which the 
fermions attract each other is a free parameter that can 
be varied continuously, and the non-interacting disper- 
sion relation is that of free fermions with an effective 
mass to*. 

The delta-shell potential can also be considered an ap- 
proximation to any central potential that is attractive 
only within a range of distances centred at r , of width 
r c <C ro, since Eq. (^) is equivalent to performing, in the 
general expression 

/>oo 

i^(|k|,|k'|) = (-iy / dr4«r 2 j l (|k|r)^(r)i«(|k / |r) J 
Jo 

(6) 

valid for any central potential (ji (x) denotes a spherical 
Bessel function), the approximation 

/>oo 

/ dr47rr 2 j,(|k|r) V (r) j,(|k'|r) 
Jo 

« r c 4^r 2 ji (|k| r ) V (r„) ji (|k'| r ) (7) 

which corresponds to taking the limit r c — > while keep- 
ing V (ro) x r c = constant = —g (g thus has dimensions 
of energy x length). A particularly simple example of 
this is the square well of Fig. |l|. 

The two-body problem associated with the delta-shell 
potential is very well-known (see Refs. [lfjjlj], for ex- 
ample). In particular, it can bind a pair in free space 
with any value of I — 0,1,2,... To simplify mat- 
ters, we will assume that the attraction takes place 
between particles with opposite spins. Finally, in k- 
space the delta-shell potential is given by V (k — k') = 
—gAirrQ sin (|k — k'| ro) / |k — k'| ro, from which it is evi- 
dent that it reduces to the contact potential in the limit 
ro — > (keeping g47rr§ equal to — Kq). Interestingly, the 
delta-shell potential, for any finite ro, does not display 
the ultraviolet divergences affecting the contact poten- 
tial. 

We will study the BCS to Bose crossover in this new 
model using the standard functional-integral technique 
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of Refs. |2 1| , |2 2| , |30|j3 1| . At zero temperature, it implies a 
description of the system in the saddle-point appropsima- 
tion, which amounts to using the BCS ground stated (as 
in Refs. |,|,|o|Jl|,|l|, for example). Thus our results for 
the ground state will be approximate, but of variational 
significance. At the critical temperature, Gaussian fluc- 
tutations about the saddle point are.-taken into account 
as in Refs. [l5|,^l],[22] (as is well-knownE3 this is equivalent 
to the RPA-like diagrammatic technique introduced by 
Nozieres and Schmitt-RinkO) . This approach is rather 
limited in that, in the strong-coupling limit, it neglects 
interactions between the preformed pairs, existing above 
T c , and so it can only describe the effect of fluctue 



on the superconducting instability at low densitiesEBEa. 
Nevertheless as we shall see it is enough to discuss the ro- 
tational symmetry breaking in the weak-coupling limit, 
in which the fluctuations are negligible, as well as the 
mechanism by which the critical temperature becomes 
larger for s-wave pairs in the BE limit (at strong cou- 
pling and low densities) . The application of these stan- 
dard methods to the DSM is fairly straight-forward, so 
we will quote here only the key expressions; details can 
be found in Ref. ||. 

Our model has four parameters: the distance tq at 
which the attraction takes place, the coupling constant 
g, the density of fermions n and their mass to*. In 
principle, the BCS ground state can be characterized 
by the dependence of the zero-temperature gap func- 
tion Ak and chemical potential fi on these four param- 
eters. Likewise, the superconducting instability can be 
described by giving the critical values of the tempera- 
ture T c and chemical potential /i c in terms of ro, g, n, m*. 
However the DSM has the remarkable property that the 
four energies Ak, /i, fcsT c ,/i c , rescaled by the "locali- 
sation energy" So = h 2 /2m*rQ (which we will denote 
A k 7 A, T c ,jl c , where k = r^k), are functions of only 
two parameters: the "dimensionless coupling constant" 
g = (eo^o) x g and the number of fermions in a sphere 
of radius tq: h = (47r/3) Tq x n. 



III. GROUND STATE 

The kernel K v (|k| , |k'|) in Eq. (0) is given by 

K l (\k\,\k'\)=-g47rr 2 (-l) l j l (\k\)j l (\k'\) (8) 

Thus although evidently the delta-shell potential is cen- 
tral and therefore not separable in the sense of the NSR 
potential, it can be written as a sum of separable terms, 
each one corresponding to a different value of I. Accord- 
ingly the gap function A k , at the saddle point, has the 
following form: 



1=0 m=-l 

In terms of the amplitudes Aj )fn the usual "gap" and 
"density" equations, for a homogeneous, stationary, non- 



polarised state with singlet pairing (A/ m = for odd I), 
read 




gMm,l'rn' fkj 



{2-kY 2£ k 



6vr2 { e c 




where e k = k 2 - fi, £ k = J e k 



A t 



and 



(4tt) 2 ji 



3V 



(|k|) y* n (k) Y Vm , (k) 



To assess the relative stability of different solutions to 
the self-consistency problem (p|JlO|), corresponding to the 
same values of g and h, one has to compare the corre- 
sponding ground-state energies: 



U = 



6n 2 



* -E 



A 



l.rn 



Ylixg 



3„. 2 



(11) 



Evidently Eq. (||) is an infinite system of non-linear 
integral equations with, presumably, an infinite number 
of non-trivial solutions and there is no systematic way 
of finding all of them. Nevertheless a certain subset, se- 
lected by the requirement that all but a few of the A; m 
are equal to zero, can be explored systematically. Since 
the effective coupling constant in the weak-coupling limit, 
Eq. (0), is (for even I) 



Ki 



-gAirrlji (fc_pr ) 



(12) 



it is clear that for small g, and within the range of den- 
sities for which Uftq < 5, we can restrict our attention 
to the first two values of I — 0,2 (see Fig. ||). These 
are the two-body stat es . wi th lowest energy (existing at 
g > 2,10, respectivelyt30) , and therefore this simpli- 
fication is also valid for our purposes in the BE limit. 
Moreover, for simplicity we will consider <i-wave states 
with a particular symmetry, choosing d^—y^iijjjeh-has 
been extensively studied in other modelsQll3EilE3llJil^-Qn 
account of its relevance to cuprate superconductivityEa 
For our two trial ground states the gap function has 
the following form, respectively: 



A s (k) = A s j Q (|k|) Y f 



(1*1) 



(13) 
(14) 



^/2 



Y 2 , 



Y 2 ,-2 (k) 



Note that this is a more complicated k-dependence than 
that of gap functions arising from interactions of the form 
(||), which depend only on the angle k. In particular, the 
gap function can change sign as k increases in the radial 
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Figure 2: The strength of the attraction in the BCS limit for 
pairing with the first four even values of the angular momen- 
tum quantum number, I = 0, 2, 4, 6. 



direction of increasing |k|, not just as the angle k is var- 
ied: see Fig. [| The oscillatory behaviour as a function of 
|k| can be regarded as a direct consequence of the singling 
out of a particular distance by the attractive interaction 
(H). More generally, we expect these oscillations, of fre- 
quency ~ l/^o, to be a generic feature of interactions 
that are attractive predominantly at some finite distance 

r °- j— i i-i 

Substitution of Eq. (pj) (Eq. (fjj)) into the self- 
consistency problem (|s|,|10D yields a much simpler prob- 
lem, which can be solved numerically for every value of 
g and n. 

For low values of g, we find that p 3> A s {fi ^ 
Kd x2 y2 ). This is the usual weak-coupling conditioncl 
characterising the BCS limit, and consequently the nu^ 
merical results display the usual generalised BCS lawo 
A s oc exp{l/NK } ( A d x 2_ y2 oc exp{l/NK 2 }) where TV 
is the free-fermion density of states, per spin, per unit 
volume. To illustrate this by an example, Fig. [| shows 
Ad x2 _ 2 vs g for constant n = 7.5. In this regime, 
the non-monotonic dependence of the effective weak- 
coupling constant Ki on the rescaled Fermi vector kpr$ 
(given in Eq. (|l^); see also Fig. ||) leads to the simi- 
larly non-monotonic dependence of A s and A x i_ y 2 on 

n = (4/9vr) (fc F r ) 3 shown in Fig. |.@ 

On the other hand, for large values of g we obtain 
p <C — A s (p -C — Ad 2 _ 2 ), which is the opposite strong- 
coupling condition, corresponding to the BE limiti. Thus 
as g is increased, while keeping n constant, p goes from 
being approximately independent of g, and equal to £p 



(= the Fermi energy ep 
having the behaviour 



k F /2m* divided by Eq), to 



(15) 

where e\ is the binding energy of the two-body bound 
state with angular momentum quantum number I = 0, 2 
(given analytically in Refs. P|[l7| for example), divided 



by £o- This evolution of the chemical potential is repre- 
sented in Fig. ||, for the d x 2_ y 2 trial ground state (the 
positive offset of p above e l b /2 that can be seen in the 
graph becomes very small, compared to e l b /2, only at 
larger values of g than those shown; additionally, it tends 
to zero as n — > 0). ~~ 

As is well knowntdtl the qualitative change of the 
ground state from BCS-like to BE-like behaviour occurs 
when the chemical potential goes below the bottom of the 
band i.e. p = 0. Fig. shows two superimposed "charts" 
of the crossover, for the s and d x i_ y i ground states, ob- 
tained using this condition. The charts include also two 
additional boundaries for each trial ground state, corre- 
sponding to p = A s ,Ad x . 2 y2 and p = -A s ,—Ad x2 _ y2 , 
which indicate the extent of the "crossover region" be- 
tween the BCS and BE limits. These charts are very 
similar to the ones presented in Ref. [ll] for s-wave pair- 
ing via the the NSR and Gaussian potentials, suggesting 
that the density-driven crossover behaviour described in 
that reference is generic to continuum models. The main 
difference that we observe for £^2^2 -wave pairing is the 
enlarged BCS region at low densities, due to the higher 
value of the coupling constant required for a two-body 
bound state. The oscillations of the boundary between 
the BCS and crossover regions at high densities are di- 
rectly related to the non-monotonic densitiy-dependence 
shown in Fig. @. 

Evidently, Fig. || suggests that at intermediate den- 
sities, at which A^ 2 _ 2 ^> A s , the energy of the trial 
ground state with I = 2 is lower than for I — 0. The 
precise value of the density at which this breaking of the 
rotational symmetry takes place, and the higher value at 
which the symmetry is restored, are given, in the limit of 
small g, by the first two positive solutions of the following 
equation: 



jo (kp) = h 



(16) 



where kp = kpr . These can be determined from Fig. |2[ 
On the other hand, for large g the system is always in 
the BE regime, in which the energy (|l l] ) takes the form 



,-/ 3 -2 



(17) 



Since the Hartree term — 3/4gn 2 is independent of I, at 
first sight this equation suggestSj-that the I = trial 
ground state, for which e\ is lowers, must have lower en- 
ergy, however note that in general lim iU s — Ud 2 _ 2 ) ^ 



YanUs — limi/d^ 2 (where U s and Ud^ ^ are the ener- 
gies of the two trial ground states and the limit refers to 
taking p -C — A s , -A^ 2 _ j2 in Eq. (|i"T|)). In fact there is 
an additioanl positive contribution to the energy, similar 
to the positive offset of the chemical potential, with re- 
spect to £j/2, seen in Fig. [| which does not appear in 
Eq. (^) because it varies slowly with g and therefore be- 
comes negligible for sufficiently large g (just like the offset 
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Figure 3: The zeroes (dashed lines) and sign ("+" and "-" symbols) of the gap function Ak on the (k x ,k y ) and (k x ,k z ) planes, 
for (a) the trial ground state with s symmetry and (b) the one with d x 2_ y 2 symmetry. 
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Figure 4: Evolution of the amplitude of the gap function for 
the d x 2_ y 2 trial ground state, as a function of g, for fixed 
n = 7.5. 



Figure 5: Evolution of the amplitude of the gap function for 
the s and d x 2_ y 2 trial ground states with increasing value of 
n, for fixed g = 2. 



of jl) . This repulsion is different for pairs with different 
internal structures, and so it is only in the n — > limit 
which Eq. ( |r7| ) allows us to conclude that the I = state 
is preferred at high g. At finite densities, the energies 
have to be evaluated numerically. Nevertheless the result, 



shown in Fig. M, confirms our expectations: between the 
two densities given by Eq. (16) the d x 2_ y 2 trial ground 
state is more stable, thus breaking the rotational sym- 
metry of the system, but only for relatively small values 
of the coupling constant. As g is made larger, the range 
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Figure 6: Evolution of the chemical potential for the d x 2_ y 2 
trial ground state, as a function of g, for fixed n = 7.5. 
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Figure 7: "Chart" of the BCS to Bose crossover for the trial 
ground state with s pairing (solid lines) and the one with 
d x 2_ y 2 pairing (dashed lines). The thicker lines are where 
the chemical potential goes below the bottom of the band, 
while the thinner lines give an indication of the extent of the 
crossover region (see text). 



of densities in which this symmetry is broken becomes 
progressively smaller until, above some critical value of 
g ~ 14, the system prefers the s state at all densities. 
Thus the region in parameter space in which the rota- 
tional symmetry is broken is relatively small. In particu- 
lar, it is confined to the BCS side of the crossover diagram 
i.e. p, > everywhwere inside the d x 2_ y 2 region.E^I 




Figure 8: Phase diagram of the relative stability of trial 
ground state with s and d x2 _ y 2 pairing symmetry. The dashed 
lines indicate the position of the phase boundary in the <? ^ 
limit, given by Eq. (|l6|). 



integral representation of the partition function 
We shall now proceed following this approach and keep- 
ing only the lowest significant corrections to the mean 
field theory. Namely, we start with a Grassman path- 
integral representation of Z for the electrons, imnlffljBsnt 
the usual Hubbard-Stratonovich transformationE3t!a£Zl to 
a functional integral over a complex pairing field A and, 
finally, expand the effective action for the fluctuations, 
Sb [A* , A] , about the saddle point of the fun|Gtk«t, 
gral above T c to quadratic (Gaussian) ordenl3li3'Eill 
This pis a well-tied approximation for the problem at 
handE2l and therefore suitable for studying the effects of 
pairing fluctuations on T c in our particular model. In 
short, using Eqs. (QJ|) to write the Hamiltonian as 



H-fiN = 



k,<r 



Z,m,q 



(where L 6 is the sample volume and Vi = 
(— 1) +1 g (47rro) 2 ) suggests that we introduce bosonic 
fields Aj jm) q (uj y ) conjugate to the operators 



Z,m,q 



5> (|k|r )y IlTn (k) c+ 2+kT c+ 2 _ ki (19) 

k 

^2 3i (l k |r )y;* m (k) c q/2 _ ki c q/2+kT (20) 



IV. CRITICAL TEMPERATURE 

Unlike the theory of the ground state that of the equi- 
librium phase at finite temperatures does not follow from 
the usual BCS theory when the superconducting insta- 
bility corresponds to the BE condensation of "preformed 
pairs" (PP). To describe such situation one must go 
beyond the mean-field theory and include fluctuations. 
This is most readily done within the framework of a path 



which evidently create and annihilate, respectively, a pair 
with total momentum q and angular momentum quan- 
tum numbers l,m. As usual, the momentum and fre- 
quency dependence of the fields captures the dynamics 
of the bosonic degrees of freedom. The additional I, m- 
dependence reflects the fact that our explicit interaction 
potential can bind pairs with different internal structures. 
Obviously the A; m of the previous section correspond 
to a homogeneous, stationary configuration of the fields, 

A; >m] q (U!„) = <5q,oAi im . 
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Proceeding in the usual wa; 



we obtain 



S b [A*,A] = E A U q K)x 

qv l.V m,m' 
xT LL,l>,rn> (q.^Af.m'flk) (21) 

where the sum on 1,1' extends only over values of the an- 
gular momentum quantum number with the same parity 
(both even or both odd), (3 = 1/fcsT is the inverse tem- 
perature and the u) v = 2vir/f3 are bosonic Matsubara 
frequencies. To further simplify the problem, and facil- 
itate the discussion of the BE limit, we follow, the pro- 
cedure employed by Zwerger and coworkeraiaEil to write 
a low-frequency, low-momentum expansion of the inverse 
propagator for the preformed pairs: 

rri r , m < (q, wv) / [di 08, m) 08, m)] 1/2 - 



E 



2 ^2 



2m, 



E 



i—x,y,z l,m,l,m 



(l-<Ji,r*m,m') (22) 



Thus after appropriate rescaling of the bosonic fields the 
known functions m ; '^ v m , fi) and /if (/3, //) play the 
role of effective boson masses and chemical potentials, as 
in Refs. ^, ^5p| , gg , for example. These two functions, 
and the rescaling factor di(/3,fj), are given in the ap- 
pendix. 

Note the different chemical potentials for bosons with 
different values of the angular momentum. Moreover, the 
anisotropic dispersion relation given by v m , {J3, fi) 
can describe not only the "rigid" propagation of a boson 
without changing its internal state, but also changes in 
its internal angular momentum through the off-diagonal 
terms, with l,m ^ l',m'. However in the BE limit, 
which as usual corresponds to [i/3 — > — oo, we have 



b.ij 

1 

l,m : l' ,m 



for I, 



and there- 



fore in what follows we shall ignore these off-diagonal 
terms (for l,m = l',m', on Jjhfi^jtilfer hand, we re- 
cover the expected behaviourEloEiH 



m Lm (A M) 



m 



b,ii 

Lm.Lm 



(AM) — 2m*). 



As usual, T c is determined by the BE condensation 
condition /if (/3, /t) = 0. This gives a different critical 
temperature for each value of I. On the other Jhand T c 
is degenerate in m = — I. . . I, as in BCS theoryH. Since 
the present method can only describe an instability of the 
normal state, our philosophy will be to compute the I = 
and I = 2 critical temperatures and then take the high- 
est of the two as the critical temperature of the system. 
Moreover we will assume that, near the critical temper- 
ature, only fluctuations with the appropriate value of I 
have to be taken into account. This is only adequate if 
the I = and I = 2 critical temperatures differ consid- 
erably, which as we shall see is the case in the BE limit. 



Of course in the opposite, BCS limit the fluctuations can 
be neglected completely. 

Under the above assumptions the "T c " and "density" 
equations take the form 



1 

9 



(-1)' 



1 



n 



Ji 



E^> 



2/ (> k ) 

2*£ 



(23) 



(24) 



where no 



3tt 



Jo 



f f/^kj i s the density of 



fermions that are unpaired above T c and the additional 
contribution coming from Gaussian fluctuations is made 
up of terms of the form 



3t jo 



rf|q||q| 2 3 /3^-U(|q|) (25) 



which correspond to fermions bound in PP with angular 
momentum given by l,m. The notations / (x) and g (x) 
have been used for the Fermi and Bose distributions func- 
tions, respectively. The "weight" w (|q|) is given by 



1/2 



tf(|q|) = n 



~ b,i 




(26) 



where the dimensionless functio 



n A ; (ji, p\ 



and the fac- 



tor k\ are defined in the appendix and each "primed" 
represents differentiation with respect to ji. This weight 
becomes unity in the BE limit /t/3 — > — oo, corresponding 
to bosons of mass 2m* each. 

Numerical solution of the self-consistency equations 
( p3| , p4] ) for f3 and ji, at the relati\jelji|rio»L| value of the 
density n = 0.5, shows the expectedlaOEiO smooth evo- 
lution between the BCS and BE limits, analogous to the 
one seen in the ground state: see Fig. ||. In particular, we 
find that the critical temperature for angular momentum 
quantum number I is 



f\ <xexp{l/7ViQ} 



(27) 



for small g but saturates to a constant value given by 



k B T c w 3.315 



2m* 



2(2Z + 1) 



2/3 



(28) 



in the large-cji limit. This asymptotic behaviour follows 
quite generally from the self-consistency Eqs. (^,^J). 

Notably, Eq. ( |2g| ) differs from the standard result (Q) 
in the presence of the degeneracy factor 1/ (21 + 1) mul- 
tiplying the density of bosons n/2. For an instability to 
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Figure 9: The critical temperature for an instability to a su- 
perconducting state with I = (solid line) and I = 2 (dashed 
line), as a function of g, for fixed n = 0.5. The increasing 
dotted lines are obtained by neglecting the contribution of 
Gaussian fluctuations to the total density, i.e. the second 
term on the right-hand side of Eq. jii]), while the constant 
dotted lines are the BE condensation temperature given in 
Eq. 
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Figure 10: The critical temperature for an instability to a su- 
perconducting state with / = (solid line) and I = 2 (dashed 
line), as a function of n, for fixed g = 3.5. The dotted 
lines rep resent the BE condensation temperatures given by 
Eq. Q. 



an ,s-wave superconducting state, with I — 0, Eq. (|2S 
reduces to Eq. (Q) and thus our result for the DSM coin- 
cides-iwith those obtaiaed earlier for models featuring the 
NSRB and contactcilO'Ea potentials. On the other hand, 
the degeneracy of the I = 2 bound state means that, at 
T c , five Bose gases condense simultaneously, but inde- 
pendently, leading to a much lower critical temperature. 
This is in contrast with the result for the anisotropic po- 
tential of Ref. [l5|. On the basis of this we conclude that 
the I = state always has higher critical temperature in 
the BE limit. 

On the other hand, for small values of g (~ 3.5) we 
find a non-monotonic density-dependence of the s- and 
d-wave critical temperatures similar to the one that we 



Figure 11: The fraction of fermions that are bound into PP 
just above T c for an instability to a superconducting state 
with I = (solid line) and I = 2 (dashed line), as a function 
of g, at the density of Fig. fil 




Figure 12: The effective mass of the PP existing just above T c 
for an instability to a superconducting state with / = (solid 
line) and I = 2 (dashed line), as a function of g, for the density 
of Fig. []. For the case of an instability to a superconducting 
state with I = 2, only the heaviest and the lightest of the 



masses m ; '* n (corresponding to 
respectively) have been plotted. 



and |m = 2 and 0, 



described for the amplitude of the gap function in the 
respective trial ground states: see Fig. n]. In particular, 
note that there is an intermediate range of densities for 
which the rf-waxe_ critical temperature is the highest. 

As expectedHlla-cil, the evolution from the BCS to the 
Bose limits is also evidenced in the fraction of fermions 
that are bound in PP just above T c , , 5n/n (with Sh = 
^ m &i jm ), and in the effective mass of such PP, 
(see Figs, [ll] and |l2|, respectively): both are negligible 
for small g, while in the large-g limit we have 5n ss n and 

m l,m ~ 2m • 

Like any theory based on a Gaussian expansion, the 
present one displays a non-monotonic behaviour of the 
critical temperature as a function of g in tl^rintermedi- 
ate regime (see Fig. W. Such enhancemenffl'E2rE2l is not 
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Figure 13: The critical value of the chemical potential for 
an instability to a superconducting state with I = (s-wave, 
solid line) and I = 2 (d-wave, dashed line), as a function of g, 
for fixed n = 0.5. The dotted lines indicate the Fermi energy 
ef and the two-body binding energies per particle e l h /2. 



presentEi in the self-consistent theoryEl due to Hauss- 
mann, suggesting that it is an artifact. It can be un- 
derstood in terms of Eq. ( p8|) and Fig. |l^ as a result of 
the PP getting lighter (m,'^ ^ 2m*) as the value of g 
is reduced. In the self-consistent theory, at least for a 
model based on the contact potential, repulsive interac- 
tions between the PP overcompensate for this, leading to 
a monotonic dependence of T c on the coupling constantEil 
(for a more advanced treatment of these repulsive forces 
see Ref. |38|; alternative methods to describe phase fluctu- 
ations and strong pairing correlations in superconductors 
are described in Refs. |3^ and Such interactions are 
completely neglected in thepresent treatment as is evi- 
denced for example in Fig. [13] which shows the evolution 
of the chemical potential, lacking a positive offset like 
the one we found in the ground state (compare the I = 2 
curve of Fig. |l3| to Fig. g). Moreover, the description of 
the Gaussian Actuations afforded by Eq. ( |22| ) turns out 
to be valid only for densities below the first maximum of 
T c , as a function of n. At the maximum, the mass of the 
PP (given in Eq. (A9), below) becomes negative, thus 
making Sh diverge. The present treatment is therefore 
only valid at small values of the coupling constant, for 
which the fluctuations can be neglected (as in Fig. [l(]: 
the only part of the plot that shows a significant con- 
tribution from fluctuations is at densities well below the 
first maximum of T c ) or at low densities, which are be- 
low the first maximum for all sizeable values of g (as in 
Fig. |). 



V. CONCLUSIONS 

We have studied exotic pairing in the context of a sim- 
ple model featuring fermions in a continuum with an ex- 
plicit, non-retarded, central interaction potential V (r): 
the Delta-Shell Model (DSM). Its novel feature is that 



the interaction is attractive only at some finite distance 
r . Because of this it provides, to our knowledge, the first 
explicit example of BCS pairing with angular momentum 
quantum number I > via a spherically symmetric (cen- 
tral) attraction. By using a standard functional integral 
approach, we have studied this breaking of the rotational 
symmetry in relation to the BCS to Bose crossover. 

By considering two trial ground states, with s and 
d x 2_ y 2 symmetries, we have found that the ground state 
with broken rotational symmetry is separated from the 
BE regime by a quantum phase transition, in which the 
symmetry of the superconducting order parameter is in- 
creased. This is due to the higher energy of two-body 
bound states with I > 0, and so it can be extrapo- 
lated to any central interaction potential. More gener- 
ally, for other models (such as those in which the single- 
particle dispersion relation and the interaction potential 
are anisotropic), our analysis suggests that a two-body 
ground state with d x 2_ y 2 symmetry is required in order 
for pairing to take that form in the BE regime. Such 
scenario is realised, for example, in a lattice model with 
nearest neighbour (n.n.) attraction and large next near- 
est neighbour (n.n.n.) hoppingo. 

Similarly, the critical temperature for superconductiv- 
ity with angular momentum quantum number I = is 
found to be higher than for I = 2 in the BE limit (of 
strong coupling and low densities). However, interest- 
ingly, this is due not to the higher energy, but to the 
related higher degeneracy of the two-body bound state 
with 1 = 2. Thus together these two observations place 
severe constraints on any interaction potential V (r) lead- 
ing to pairing with I > in the BE limit. 

In our model, the rotational symmetry-breaking is a di- 
rect consequence of a non-monotonic dependence of the 
superconducting properties on the fermion density which 
is present only in the BCS regime. Such rise-and-falls 
can be understood in terms of the oscillatory form of 
the "gap function" in k-space, whose frequency is ~ Tq , 
and presumably they are generic to interaction poten- 
tials that are attractive predominantly at a finite dis- 
tance. In Refs. []l]j4l] the possible implications of our 
model to cuprate superconductors, on the basis of the 
similar behaviour observed in the doping-dependence of 
the supercondcuting gap and the critical temperature, 
were discussed. A similar rise-and-fall has been known 
for some time in nuclei (see Ref. for example). Of 
particular interest, in connection with recent theoreti- 
cal speculations on supepfluiditw-jn magnetically trapped 
gases of fermionic atomsc3oc3o (for an informal review 
and further references see Ref. [I?]) , is the possibility that 
the present mechanism would lead to exotic pairing for 
sufficiently high densities in these systems. Interestingly, 
because the change from s- to d-wave pairing is a quan- 
tum phase transition, it can take place at arbitrarily low 
temperatures. On the other hand, the density would have 
to be raised until r s ~ r (where, in order to achieve a 
phenomenological description of the inter-atomic poten- 
tial, ro may be taken to be roughly the size of a diatomic 
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molecule) . 

In relation to possible future work, we end by not- 
ing that we have described the ground state of the DSM 
in terms of a homogeneous saddle point, and we have 
only taken into account pairing fluctuations around that 
saddle point. In principle, by performing a more gen- 
eral Hubbard-Stratonovich transformation, including ad- 
ditional fields associated with the density (not just the 
pairing amplitude), one could study the effect of den- 
sity fluctuations as well as the possibility of phase sepa- 
ration through a first-order gas-liquid phase transition: 
physically, one expects that the attraction at a finite 
distance could favour, in addition to the pairing with 
I > which we have considered here, the formation of 
clusters of more than two particles Xas in lattice mod- 
els with nearest-neighbour attractionlla, and unlike those 
with on-site attraction^]; for a discussion of the similar 
phenomenon of "quartetting" see Ref. see also foot- 
note 19 of Ref. |3§). Evidently, this would be very inter- 
esting in the light of recent ^discussions of inhomogeneity 
in cuprate superconductorsE2l. 



Acknowledgments 

We thank K. Capelle for detailed comments on an 
earlier manuscript and A.S. Alexandrov, N. Andrenacci, 
N.W. Ashcroft, J.M.F. Gunn, P. Howell and H.M. Klein- 
ert for useful discussions. JQ acknowledges finan- 
cial support from the TMR programme (EU; contract 
No. ERBFMBICT983194), FAPESP (Brazil; process 
No. 01/1046 1-8) and CNPq (Brazil). 

Appendix A: EXPANSION OF THE INVERSE 
BOSONIC PROPAGATOR 

The inverse bosonic propagator in Eq. ( pi] ) is given by 



SlJ'Sm.m' ~ S i l l M,m,l',m> ( k ) X 

xG ( § + k, iu n j G (— - k, iw v - iuj n ) (Al) 



where Go (k, iu n ) = (iuj n — £k)~ (with £k = eqe^) is the 
free fermion Green's function and the oj n = (2n + 1) tt / j3 
are fermionic Matsubara frequencies. The derivation, 
starting from Eq. jig)), is entirely analogous to that 
of the similar expression in Ref. glj for example. In 
Ref. the full frequency dependence of r _1 (q, iu> u ) 
was taken into account to obtain the critical tempera- 
ture of a model featuring the-rx>ntact potential. The pro- 
cedure that we follow hereEiEil yields the same results 
in the BCS and BE limits and a much simpler numer- 
ical problem in the crossover regime (where any theory 
based on a Gaussian expansion must be regarded as an 
interpolation scheme anyway) . First we analytically con- 
tinue the second Green's function on the right-hand side 
of Eq. ( |aT| ) with respect to the bosonic Matsubara fre- 
quency, G (§ - k, iu v - iw n ) — > G (f - k, w - icj n ) , 
and then we perform the summation over n. Using the 
contour G in Fig. [l4| (which we deform into Gi and G2) 
we obtain 



L 3 1 

§l,l>6 m ,m> - /J 1 ~ l &l,m,l',m' ( k ) X 

q z — ' 

1 - / Cggq/2+k) ~ f(P jgq/2-k ~ H ) 
£ q/2+k + ^q/2-k ~ w 



It is now easy to write a low-frequency, low-momentum 
expansion of the form 



T l.rn,V,m> iu » ) ~ OJ.m.l'.m' (A fl) ~ idl, m ,l>,m> (fi, H) + ^ ^^ c lm,l', m ' (A M) QiQj ( A 3) 

i,j=x,y,z 



by— .simply differentiating with respect to q and 
wE5 We find that the coefficients ai im ,l',m' and 
dLm,i',m' are diagonal in l,m, and degenerate 
in m: ai, m ,i',m' (fl, /•*) = oiW^l/W and 
di, m ,i',m' (/3,m) = di(J3,fj,)di t i>6 mtm t. The dimen- 

sionless functions o; fl\ = r Q e L^ 3 ai (/3, /Lt) and 



{j3, fij = r Q elL 3 d t (/3, fj,) are given by 



ai 



(Am) = I -(-if- f 

v ' 9 1" Jo 

Xjl ( 



2 -Td 










k 




k 


71" JO 









1 - 2 



2er 



(A4) 
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Figure 14: The contours used to perform the summation over 
the fermionic Matsubara frequencies Lu n in Eq. (Al). 



(AA) = (-1)' 



x (A5) 



1 - 2 / (fit) 0f (fei) 



4e? 

k 



where we have used the notation /' (a;) = df (x) jdx and 
the dimensionless inverse temperature is (3 = £q/3. The 
coefficient in q 2 has the form 

4fm,l', m ' (Am) = \dl (A M) ^,!'5m,m'*i,j 

V CM (A6) 



where «; 



l.m.l' ,m' 



h ( l) d2 ^ Y Cm(f) *Ml>,m> (k), 



which can be evaluated easily by writing it in terms of the 
Gaunt coefficients of relativistic quantum mechanics (see 
Ref. pU for example; note in particular that k)'L > „ = 



for all i ^ j) and, finally, A/^ 
given by 



(aa) 



(a a) 



l,m,l,m 

n)£oAz,i' (/3,/i) is 



J i 



(l k l) 

^ 2 r (/%) 



2% 



(A7) 



The integrands of Eqs. (A4), (A5) and (A7) have no 
poles on the domain of integration, and therefore are 
straight- forw ard to evaluate numerically. From compar- 
ison of (A3) to ( p2| ) it is evident that the masses and 
chemical potentials of the PP are given by 



At (A j") = -h (aa)M (A m) 



(a a) 



1/2 



(A8) 
-l 

, (A9) 

where we have introduced the definitions 

A? (A AJ = eoVMAM) and TO^, m - (/3,A) = 

(2m*) nrf'mi' m i (A A*) an d we have ommitted the 
dependence of some of the functions defined above on 
/3, /i for brevity. 
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For fermions in a continuum with a pspulsive interaction 
V (r), the Kohn-Luttinger meehanismuJ can lead to pairing 
with high angular momentumlSa. But note that in this case 
the BE limit can never be realised, since obviously there is 
no two-body bound state. 

Note the similarity of this behaviour to that of some lattice 
models (c.f., for example, the dependence of the critical 
temperature on band filling shown in Fig. 18 of Ref. |^). 
Note that the fact that the phase transition at these quan- 
tum critical points is of first-order is due to our choice 
of trial ground states, Eqs. (|l^,[l4|). Had we allowed for 
mixing of the s and d x i_ y i order parameters, we would 
presumably have found two second-order (instead of one 
first-order) phase transitions, in analogoy with the similar 
results in Ref. |l2]. 

Note that Eq. ( IA2) does not coincide with the similar ex- 
pression in Ref. |22[ obtained by analytical continuation 
of the inverse bosonic propagator: rj~ m p m > (q, *w„ ) — > 

r~~m c m' (l> w )- Apart from the model-specific features, 
which are our main concern here, it differs also in the pres- 
ence of the continuous variable w in the argument of one of 
the Fermi distribution functions. At the bosonic Matsubara 
frequencies w = iov, we can write / f/3 [e q /2-k ~ ituj\ ) = 
/ (/3e q /2-k) and therefore both expressions are identical 
when the full frequency-dependence of ; , m , (q , io y) is 
taken into account. However the expression in Ref. E2| does 



not admit a small-frequeaejy expansion of the form (A3) 
because as is well knownpo it has a branch cut along the 
real axis that crosses the imaginary axis whenever /i > 0. 
(On the other hand such expression is the correct start- 
ing point for the-iderivation of a time-dependent Ginzburg- 
Landau theorycj, which obviously involves an expansion 
along the real axis.) 



